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Pressure and non-linear quark number susceptibilities in QCD 

R. V. Gavai^ and Sourendu Gupta**. 

'^Dept. of Theoretical Physics, Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 
400005, India 

We report the first measurements of sixth order non- linear quark number susceptibilities in continuum QCD. 
This extends our earlier computation of the continuum limit 1 of non-linear quark number susceptibilities and the 
pressure in QCD for m/Tc — 0.1. For T/Tc > 1.5, no phase transitions are observable for n < SOT. Systematics 
of off-diagonal susceptibilities are reported. The equation of state in the high temperature phase of the QCD 
plasma is constructed. 



The derivatives of the pressure with respect 
to the quark chemical potentials (one for each 
flavour) — 
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are the quark number susceptibilities (QNS, for 
I > 2) and the number densities (n, for 1 = 1). 
Thus, the quark number susceptibilities measured 
at zero chemical potential (where direct lattice 
computations are possible) allow reconstruction 
of the equation of state at finite chemical poten- 
tial through the usual multivariate Taylor series 
expansion — 



AP(T,M)= ^X? 
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where AP(r, = P(r, - P{T, {0}). The 
coefficients for odd I vanish due to CP symme- 
try. The convergence of eq. ||2J) is spoilt by phase 
transitions; hence any bound on the radius of con- 
vergence limits the location of the nearest phase 
transition to the point at /i = 0. Due to finite lat- 
tice spacing ambiguities in X/gi...) a^nd the rapid 
convergence of the Taylor expansion away from 
a phase transition, the best way to compute the 
continuum equation of state (EOS) is through the 
Taylor expansion [T|. We show later that, for 
T > Tc, the equation of state |2I1I3| is dominated 
by the I = 2 term (here Tc means the transition 
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or crossover point for /i = 0). The QNS are also 
interesting in themselves 

We report measurements of the continuum 
limit of the susceptibilities upto order 6 in QCD 
using the quenched approximation for T > l.STc, 
where the difference between quenched and dy- 
namical computations is expected to be of the 
order of 5-10%. Near Tc, where the difference is 
large, we show results obtained at finite cutoff in 
dynamical two-flavour computations. 

Off-diagonal susceptibilities need the computa- 
tion of quantities such as ti A" for matrices A. 
Independent sets of Ny/n random complex Gaus- 
sian vectors, {r}, {s}, etc., are used for the noisy 
evaluation such as ti A^ = (r|A|r)(s|^|s). For 
the evaluation of 0ii, where A is anti-Hcrmitean, 
the result is negative: in the Ny —^ oo limit, the 
histogram of noisy evaluations is skew, peaking 
at zero, with a tail to the left, and vanishing on 
the right. As shown in Figure ^ fairly large Ny 
are needed to get a signal (we used Im (r|^|r) 
for noise reduction) near T^. For larger T even 
Ny = 100 docs not give a signal pp. The con- 
jugate gradient stopping criterion is not a cru- 
cial parameter in this computation. Of the off- 
diagonal susceptibilities, only Xud and Xuudd were 
found to differ from zero significantly near Tc (de- 
generate 3-flavour QNS such as Xuds and Xudss 
also vanish within errors). The temperature de- 
pendence of these quantities is shown in Figure 
[5] The peaks in these quantities close to Tc are 
crucial for building up the peaks in the pressure 
near the critical end point [2], since the diagonal 
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susceptibilities go smoothly to zero near 
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Figure 1. Normalised histogram of measurements 
of 011 with different iV„ for dynamical Nf = 2 
computations at l.OSTc on a 4 x 10'^ lattice with 
m/Tc = 0.1. 
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Figure 2. Temperature dependence of Xud (left 
panel) and Xuudd (right panel). Boxes denote 
data for quenched QCD with Nt = 14, and circles 
for dynamical Nf = 2 QCD with Nt = 4. 



The 6th order diagonal QNS, Xuuuuuu, can 
be reconstructed by standard methods from the 
sixth derivative of Z — 

^(0111111 + 14011112 + 2001113 
+4501122 + 150114 + 600123 + 6015 
+ 150222 + 15024 + 10033 + 06), (3) 



where the operators other than 0e have been 
written out in In the notation of [Jj this op- 
erator is given by 

06 = -120(6 - 1) + 360(4 - 1®2) 
-210(2 • 1©2 - 2) -60(1 ©2© 1 ©2) 
-120(3 • 1 © 3) + 30(3 • 2) + 60(1 © 2 © 3) 
+60(1 © 3 © 2) + 30(2 • 1 © 4) - 10(2 • 3) 
-15(2©4)-6(l©5) + (6) (4) 

Results for the 6th order susceptibilities are 
shown in Table The large finite lattice size 
effects shown in Tableland in the other QNS ^ 
underline the need for taking the continuum limit 
in order to get reliable estimates of the location 
of the critical end point, where the Taylor ex- 
pansion diverges due to the growth of the higher 
order susceptibilities. 

Estimates of the radius of convergence of eq. 
^ at the 4th and 6th orders are 
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fj!"^^ /T w 4.6 in the continuum limit for 1.5 < 
T/Tc < 3 [2, whereas fJ.'-^^ /T « 28 for Nt = 14 
and is possibly significantly larger in the contin- 
uum. For T > Tc, fj-* limit the phase boundary 
between the plasma and the appropriate colour- 
superconducting phase. If the phase boundary 
is further off, or if, for T < Tc, there is quark- 
hadron continuity |^, then these two numbers are 
the value of fi/T at which the two terms being 
compared give equal contribution in the Taylor 
expansion. At SPS energies, where /i/Tc — 0.45 
the 4th order term gives a 5% correction to 
AP/T'^, and the 6th order term is totally neg- 
ligible for T/Tc > 1.5. At the RHIC, where 
fi/Tc — 0.15, the leading term contributes more 
than 99% of the total. 

The relation between AP{T, {fi})/T'^ and the 
quark number density n(T, {^})/T^ is the equa- 
tion of state P]. Since n is the first derivative of 
the pressure, the EOS is completely determined 
once the susceptibilities are known. As shown in 
Figure 121 the difference between the leading term 
involving Xum and the next term, becomes im- 
portant only when n/T^ « 1, i.e., for /i « T. 
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Figure 3. The equation of state in the continuum 
hmit at finite chemical potential for T = 2Tc with 
the leading term only (dashed line), and the first 
two terms (full line) of the Taylor series. The 
band shows the l-cr error due to statistical er- 
rors in the susceptibilities. The EOS at 1.5Tc and 
3Tc is indistinguishable at this scale. The dotted 
line is the response function d{P/T'^) /d{n/T^) at 
fixed fi and T — the instability in |2j is an artifact. 



The remaining terms change the predictions by 
less than the statistical error band shown. Away 
from Tc, this justifies, post-facto, the scaling of 
the EOS obtained by reweighting at = 4 to 
give the continuum limit, as performed in 0. Fi- 
nally, Figure |21 shows that the measurement of 
susceptibilities is the simplest and most accurate 
route to the EOS in this range of T and /i which 
is important for experiments. 
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